_ " | VOLUME 80 SEPARATE No. 578 


PROCEEDINGS | 


SOCIETY 


CIVIL ENGINEERS 
DECEMBER, 1954 


INTERACTION BETWEEN WELL 
AND AQUIFER 
| by Wen-Hsiung Li, A.M. ASCE 
| HYDRAULICS DIVISION 
| {Discussion open until April 1, 1955} 
| Copyright 1954 by the American Society oF Civit ENGINEFRS 
— Printed in the United States of America 
Headquarters of the Society 
33 W. 39th St. 
New York 18, N. Y. 
PRICE $0.50 PER COPY 


\ 
AMERICAN 
SOCIETY OF 
cv 
ENGINEERS 
FOUNDED 


THIS PAPER 


--represents an effort by the Society to deliver 
technical data direct from the author to the 
reader with the greatest possible speed. Tothis 
end, it has had none of the usual editing required 
in more formal publication procedures. 


Readers are invited to submit discussion apply- 
ing to current papers. For this paper the final 
date on which a discussion should reach the 
Manager of Technical Publications appears on 
the front cover. 


Those who are planning papers or discussions 
for “Proceedings” will expedite Division and 
Committee action measurably by first studying 
“Publication Procedure for Technical Papers” 


(Proceedings — Separate No. 290). For free 
copies of this Separate—describing style, con- 


tent, and format—address the Manager, Techni- 
cal Publications, ASCE. 


Reprints from this publication may be made on 
condition that the full title of paper, name of 
author, page reference (or paper number), and 
date of publication by the Society are given. 


The Society isnot responsible for any statement 
made or opinion expressed in its publications. 


This paper was published at 1745 S. State Street, 
Ann Arbor, Mich., by the American Society of 
Civil Engineers. Editorial and General Offices 
areat 33 West Thirty-ninth Street, New York 18, 
N.Y. 


| 


INTERACTION BETWEEN WELL AND AQUIFER 


Wen-Hsiung Li,’ A.M. ASCE 


SYNOPSIS 


The piezometric head along a well is not constant as usually assumed, but 
it must decrease towards the discharging end of the well. An investigation is 
made to find the effect of this variation of piezometric head along the well on 
the flow in the aquifer. It is found that the formation loss to the discharging 
end of the well involves an additional term roughly equal to the velocity head 
at the discharging end of the well. The loss through the well screen is usual- 
ly negligible. What is usually considered to be the screen loss is actually the 
additional loss mentioned above caused by the variation cf piezometric head 
along the well. A solution is obtained for the axially sy:.. aetric flow from a 
confined uniform aquifer into a fully penetrating well. 


INTRODUCTION 


In the study of flow into wells, the assumption of a uniform distribution of 
piezometric head on the surface of the well screen is generally made. Fora 
steady flow into a well fully penetrating a confined uniform aquifer, this as- 
sumption leads to a two-dimensional radial-flow solution which has been 
widely used in practice. This well known solution is 


H - Hog H H,= (1) 


in which H is the piezometric head in the aquifer, Q is the rate of discharge 
in the well, k is the coefficient of permeability of the soil (dimensionally, 
length per unit time), L is the thickness of the aquifer, Dp is the outside di- 
ameter of the well, and r is the horizontal distance from the center of the 
well (see Fig. 1). The subscripts o and d are used to signify values at the 
outside surface of the well screen and at the discharging end of the well re- 
spectively. (The letter symbols used in this paper have been collected in 
Appendix II for convenience of reference). However, by energy or momentum 
consideration, the piezometric head h in the well must decrease towards the 
discharging end of the well. It follows that the piezometric head Hp, in the 
soil just outside the well cannot be uniformly distributed along the screen. 
The purpose of the present investigation is to find, through a theoretical an- 
alysis based on wel! established principles, the effect of this variation of 
piezometric head along the well on the flow in the aquifer. A steady flow 
from a confined uniform aquifer into a fully penetrating well is used for 
demonstration. 


1. Asst. Prof. of Civil Eng., The Johns Hopkins Univ., Baltimore, Md. 
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The Flow in the Well 


In the flow in the well, the piezometric head h, the rate of discharge Q 
and the rate of flow added to the well through the screen dQ/dz, v.-y with z 
along the well. The vertical distance z is measured from the bottom of the 
well as shown in Fig. 1. The equations governing the flow are the momen- 
tum equation, 


and the orifice equation for the screen openings, 


aa 29(H.- h) = 2g(H,-h) (4) 


in which A is the inside cross-sectional area of the well, g is the gravitation- 
al acceleration, C, is the coefficient of discharge for the screen openings, 
Ap is the screen opening area per unit area of screen surface, Cg = CcAp is 
the screen coefficient, and D is the inside diameter of the well. In these 
equations, the quantities h, Q, dQ/dz and Hy vary with z. In the momentum 
equation, the friction on the screen-wall has been neglected. The screen co- 
efficient Cs in Eq. 4 can be determined experimentally. The validity of Eqs. 
2 and 4 has been demonstrated experimentally by Petersen, Rohwer and 
Albertson.2 The principle of conservation of energy is not used here since 
the loss of energy, mainly due to the impact between the incoming water and 
the water already flowing in the well, is unknown a priori. The loss of energy 
in the flow, however, can be determined from Eq. 2. 

The flow in the aquifer must be so adjusted that these three equations are 
satisfied. Substitute into Eq. 2 the value of Q from Eq. 3 and the value of h 
from Eq. 4, 


L L 
dQy_ 2Q,/ dQ 
He hg 29(C.xD)* Ge) oh dz az] 


Eq. 5 is the condition that the values of Hp and dQ/dz must satisfy at every 
point z along the screen. 


2. Petersen, Rohwer and Albertson, “Effect of Well Screen on Flow into 
Wells”, Proc. ASCE Vol. 79, Separate No. 365, Dec. 1953. 
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The Flow in the Aquifer 
In this investigation, attention is confined to cases with laminar flow 


throughout the aquifer. For these cases, the flow throughout the aquifer fol- 
lows Darcy’s Law for seepage: 


in which qy and qz are the volume of flow per unit time per unit cross- 
sectional area of soil in the r- and z-direction respectively. For continuity 
in this axial-symmetric flow, the distribution of the piezometric head H in 
the aquifer must satisfy the following equation:3 


r or or? 0 


A solution of Eq. 7 for the steady flow in a uniform aquifer, satisfying Eq. 5 
at the screen of a fully penetrating well, can be obtained for cases satisfying 


the following condition: 
acD, 
kD, 
A mD 


This condition is generally satisfied in practical cases. The solution is (see 
Appendix I) 


in which 


g 


3. Muskat, M., “Flow of Homogeneous Fluids Through Porous Media,” 


McGraw Hill Co., 1937, p. 141. 
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4 
SC 
H= Hoa + + ight, n*(I+N,) 
= 
oo (9) 
8 (-1)" nnz 
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and, for practical purposes, 


since the contrioution of N,, to the sum of the series is negligible in practical 
cases. In these equations, "Ko and Kj are the nodified Bessel Functions of 
the second kind »f zeroth and first order respectively. Both K, and Kj de- 
crease with an increase of the argument in a manner similar to the exponen- 
tial functions. Eq. 9 gives the three-dimensional distribution of piezometric 
head in the aquifer. The flow along the top of the aquifer experiences a 
greater loss of head than that along the bottom. The difference can be found 
from Eq. 9 to be Qj/(gA”) and is shown schematically in Fig. 1. 

The distribution of flow in the aquifer can be obtained from Eq. 9 accord- 
ing to Eq. 6. The distribution of inflow to the well, dQ/dz, can be shown to 
be 


xD,k(2H) = 4 K ng Ds -cos("E) 


Since the last term in the parenthesis does not amount to more than 10° 2 in 
practical cases, 


dQ. Qs (11) 
dz a 
This shows that the velocity distribution is practically the same as the radial- 
flow approximation although the distribution of piezometric head near the well 
is quite different from the approximate solution (Eq. 1). The explanation for 
this is that, due to the high resistance offered to the flow by the aquifer, a 


slight change of velocity distribution would cause an appreciable change of 
piezometric head distribution. 


Practical Significance 


Although the actual distribution of the piezometric head or of the velocity 
of flow is of little practical interest, Eq. 9 yields interesting results with re- 
gard to the formation loss and the loss across the screen. 


Formation loss to discharging end of well 


The head loss in the aquifer from a point at distance r to the discharging 
end of the well is (H-H,q). As r increases, the last term in Eq. 9 rapidly 
diminishes. Thus for reasonably large value of r (say 2L or larger), the 
value of (H-H,,4) is 


H-H,4~ loge) + (r>2L) (12) 


2xkL 


co Qo 
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It is thus shown that the formation loss to the discharging end of the well is 
greater than that obtained from the approximate radial-flow solution (Eq. 1) 
by an amount equal to 4/3 times the velocity head of the flow in the well at 
the discharging end (see Fig. 1). 

In this analysis, only cases with laminar flow throughout the aquifer are 
considered. For cases with turbulent flow near the well, there will be an ad- 
ditional loss of head. 


Head loss across the screei 


To find the difference in piezometric head across the screen, Eq. 4 may 
be reduced to read 


H.-h = 


Substitute Eq. 11 into Eq. 13, 


H.-h = (13a) 


It can be seen from this equation that, with reasonable values for the screen 
coefficient2, the piezometric head difference across any point of the screen 
amounts to only a very smail fraction of the velocity head of the well dis- 
charge. 

Thus the total loss of piezometric head to the discharging end of the well 
is 


hy fog, (20) 4 4-5 + (2 
H hg loge| 3 * lac) (r=2L 


This is the sum of the formation loss along the top of the aquifer and the loss 
through the screen. The same loss of piezometric head can be obtained along 
any other streamline. For example, the formation loss along the hottom of 
the aquifer to the bottom of the well is less than that along the top by an 
amount equal to / (g& ).(from Eq. 9, see Fig. 1), but the loss of piezometric 
head inside the well from the bottom to the top is equal to Q3/ (gA ) (from Eq. 
2), thus making the total loss to the discharging end of the well the same as 
that of the flow along the top of the aquifer. 

When no turbulent flow in the aquifer is involved, the data obtained from 
well tests are usually expressed as 


H -h, = aQ, + (15) 


where a and b are constants for a given well and given r. The term aQ, is 
usually referred to as the formation loss while the term bQ3 is considered to 
be the screen loss. By comparing Eqs. 14 and 15, athe term aQ, does not rep- 
resent the entire formation loss, and the term bai is mainly that part of the 
loss caused by the variation of piezometric head along the well. 
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CONCLUSIONS 


Although the analysis presented in this paper is confined to the case of a 
fully penetrating well in a confined aquifer, the following conclusions are be- 
lieved to be applicable to wells in general. 

1) With laminar flow throughout the aquifer, the formation loss to the dis- 
charging end of the well consists of two parts. One part is of the same order 
of magnitude as the velocity head in the well at the discharging end. (For a 
fully penetrating well, this part is equal to 4/3 times the velocity head). The 
other part is equal to that obtained with the assumption of uniform piezomet- 
ric head at the well surface. 

2) The velocity distribution in the aquifer obtained by assuming uniform 
piezometric head at the well is accurate enough for practical purposes. 

3) Since the velocity distribution obtained by assuming uniform piezomet- 
ric head at the well is accurate enough for practical purposes, Eq. 13 may be 
used for finding the difference of piezometric head across the screen with the 
approximate value of dQ/dz thus obtained. This head loss across the screen 
is negligible for screens in reasonably fair condition. What is generally con- 
sidered to be the screen loss in Eq. 15 is actually the additional loss caused 
by the variation of piezometric head along the screen. 
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APPENDIX I - SOLUTION OF EQ. 7 FOR WELL SHOWN IN FIG. 1 


Let = kHL/Qg, x= r/L, y= 2/L, 5=Do/2L, QgkDo/(g#) and 
B = (D/CsL)*. Eq. 7 can then be rewritten as 


ay * x ax * axe 


For the aquifer and well shown in Fig. 1, the boundary conditions to be satis- 
fied are 


2 = 0 at y= 0, since qz must vanish there; 


2) 32 = 0 at y= 1, since qz must vanish there; 


3) ¢ —— as X —- co, Since at large di.‘tance from the well the flow is es- 


sentially radial; 
—0asx-—-oco , for the same reason; 
=%q at x= 6 and y = 1, by definition; and 
6) For any point along the well screen, Eq. 5 must be satisifed. 


7 
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Since by continuity of flow across the screen, 


Eq. 5 may be rewritten as 


2 


Since at y = 1, the two terms involving an integral vanish, and ?9= ¢oq and 


- od by definition, 


Thus B.C. (6) is that, for any value of y at x = 6, the following condition must 
be satisfied: 


[ (22) -(22)'] 


2 
2nd y— 2n KS 
J y 


= Fcx)- 
Eq. 16 can then be written as 


(22) 
F 


where c must be a constant. Furthermore, to satisfy B.C. (1) and B.C. (2) 
simultaneously, no solution is possible if c is an imaginary number. There- 
fore c must be a real constant. Here the primes and double-primes signify 
the first and second derivatives of the functions respectively. Two cases 
must be considered separately, with c= 0 and with c 0 respectively. 
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= nmD.k ($2) (17) 
Qe 
(20) 
Let 
- (21) 


1) When c = 0: 
When c = 0, G''= 0. Then G' = Ap where Apo is a constant. In order to 
satisfy B.C. (1), (2) or (4), Ag must be zero. Therefore 


G =B, (23) 


where Bo is a constant. With c = 0, Eq. 22 gives 
F=A,+ B, loge X 
where A, and By, are constants. Thus, with c= 0, 
FG= a+ logeX 


where ap and by are constants. 


2) When c 7 0: 
When c # 0, Eq. 22 gives 


G = A, cos(cy) + Bz sin(cy) 


G'= -cAzsincey) + ¢B, coscey) (27) 


To satisfy B.C. (1), G' must vanish at y= 0. Therefore cBy =0, or Bg = 0 
since c # 0. To satisfy B.C. (2), G' must vanish at y = 1. Therefore 

cA: sin(c) = 0. Since c ¥ 0, and with vanishing Ag the solution would be 
trivial, we must have sin(c)= 0. Therefore c= n7, where n = 1,2,3,......, and 


G =a, cos(nxy), (n=1,2,3,----) (28) 
where a, is a constant for each value of n. With c = n7, Eq. 22 also gives 


F= A,l.(nex) + (29) 


where I, and Kp are the modified Bessel Functions of zeroth order of the first 
and the second kind respectively. The value of Ip increases while the value of 
Ko decreases with an increase in the value of the argument in a manner simi- 
lar to the exponential functions. To satisfy B.C. (4), FG' must approach zero 
as x increases. Since G' does not vanish identically, F must approach zero 
as x increases. Therefore Ag must be zero. Thus, when c = n7, 


= F-G = ay cos(nxy), (n=1,2,3,-----) (30) 


The constant Bg in Eq. 29 has been absorbed into a,,. 
578-8 


(24) 
(25) 
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Since both Eq. 25 and Eq. 30 are solutions of Eq. 16 which is linear, the 
general solution of Eq. 16 is therefore 


oo 


$= a, + bo logex + K.(nrx)-cos(nry) (31) 


n=! 
The constant b, can be determined from B.C. (3): Using the identity 


—nn- K,(nrx) 


we have 


be _ nit On 33) 
Ox x 


B.C. (3) is satisfied with by = 1/(27), since K;(n7™x), the modified Bessel 
Function of first order of the second kind, vanishes at large values of x. The 
constant a, can be determined from B.C. (5): Atx=6 andy = 1, Eq. 31 and 
B.C. (5) gives 


n=} 


Therefore 


= Pog + loge(4) On Kind) 


+). a, 


n=| 


This solution of Eq. 16 satisfies the first five boundary conditions. The co- 
efficients a, must be such that B.C. (6) is also satisfied. 

Eq. 35 indicates that the solution consists of the solution for two- 
dimensional radial flow (Eq. 1) with additional perturbation terms under the 
summation signs due to the non-uniform distribution of piezometric head 
along the well screen. To linearize Eq. 20 (B.C. 6), cross-products of the 
small perturbation terms are neglected. For example, we have from Eq. 33 


the value ot 9¢ and we use in Eq. 20 
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Ll 
ax = (32) 
: 
~ 
~ 


fats) I’m dn Kiedy (26) 


In this manner, Eq. 20 may be reduced to the following form: 
oo 


an — on K.(nd) 


n=| 


(37) 


n=| 


co 
of 2 
n=! 

The values of the constants a, must be such that this equation is satisfied 
for every value of y between 0 and 1. A more managenble equation may be 
obtained by multiplying Eq. 37 by cos(m7ty)dy where m = 1,2,3... and integrat- 
ing term by term from y= 0 to y= 1. Since 


0,if m#n (38) 
0 +, if m=n 


if m=n 


2 


J cosimny) dy = 0 


(22) 
Ox ‘od 
. 
& 
and 
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for each value of m, we have 


[ + + < Kind) Om 


With different values of m, Eq. 42 gives an infinite series of simultaneous 
linear equations with the infinite number of am as 4owns. Since the value 
of a,, is small for large values of m, an approximate solution to the problem 
can be obtained by using a finite number of terms under the summation signs 
in Eq. 35. For this purpose, the values of n and m in Eq. 42 may be limited 
to a finite value, and the finite number of a,, can be determined from the re- 
sultant finite number of simultaneous equations. However, the process of de- 
termining the coefficients in this manner is very tedious and the result is only 
approximate. For the problem under consi:deration, the value of @ in practi- 
cal cases is extremely small, usually less than 16-3, Thus, as a first approx- 
imation, a,, may be computed from 


A second approximation may be obtained by first computing the values of a, 
in Eq. 42 according to Eq. 43: 


+ m + On = 


2 (44) 
n(m* 


It can be shown that, if 


K, 
(45) 


as in nearly all practical cases, Eq. 44 is reduced to Eq. 43: First the term 


SK (m6) is negligible compared with K)(m7‘6) for any value of m because 
of the following relationship among the Bessel Functions: 


< 2 
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: 
(42) 


Thus 


(15) 


Secondly, for the same reason, the term involving the summation sign in Eq. 
44 is negligible compared with unity: 


Ky(nnd) 


(48) 
ied) ay 
nem 


Since the condition specified by Eq. 45 is generally satisfied in practice, Eq. 
43 is accurate enough for practical purposes and is presented with Eq. 35 as 
Eq. 9 in the text of this paper. 


APPENDIX II - NOTATION 


The following symbols, adopted for. use in this paper and for the guidance 


of discussers, conform essentially with American Standard Letter Symbols 
for Hydraulics (ASA Z10.2 - 1942), prepared by a Committee of the American 
Standards Association, with Society representation and approved by the Asso- 
ciation in 1942: 


A cross sectional area of well; 

A area of screen openings per unit area of screen surface; 
Ao, Al, 49 Bo, Bi, Do, etc. = constants of integration; 

2m, an = constants in solution of Eq. 16; 

a, b = constants in well loss formula, Eq. 15; 

C. = coefficient of discharge of screen openings; 

Cg =C,Ap, the screen coefficient; 

c = constant; 

D_ = inside diameter of well; 

Do = outside diameter of well; 

d = subscript signifying values at discharging end of well; 

F, G= functions; 

g =gravitational acceleration; 

H =piezometric head in aquifer; 

h =piezometric head in well; 

I, =modified Bessel Function of zeroth order of the first kind; 
Ky =modified Bessel Function of zeroth order of the second kind; 
Kj =modified Bessel Function of first order of the second kind; 

k =coefficient of permeability of soil (dimensionally length per unit time); 
L_ =thickness of aquifer; 


578-12 


| 
m 44 
3 
1 
< 
« 


m,n = positive integers; 

Nn function defined in Eq. 9a; 

oO subscript signifying values at outside surface of screen; 
rate of discharge passing any section of well; 

discharge per unit cross-sectional area of soil; 
horizontal distance from center of well; 

r/L; 

z/L; 

vertical distance from bottom of well; 


s 
= D,/(2L); 
= kHL/Qq 


Q 
q 
r 
x 
y 
a 
B 
6 
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= QykDo/(gA*); 
D 2 : 
= Gz) 
| 
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